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Abstract 


A  simple  model  is  developed  which  permits  us  to  calculate  the  sudden 
ionospheric  frequency  deviation  which  would  be  caused  by  a  known  solar 
X-ray  burst,  or  vice-versa,  to  calculate  the  time  profile  of  the  X-ray 
burst  responsible  for  an  observed  frequency  deviation.  This  model  is 
also  used  to  compare  the  frequency  deviations  observed  during  six  solar 
flares  with  the  time  profiles  of  the  solar  emissions  at  centimeter  and 
X-ray  (0.5-10  k)  wavelengths  which  accompanied  the  flares. 

Better  time  and  spectral  resolution  are  needed  jn  order  to  permit  a 
detailed  comparison  of  the  variations  in  the  X-ray  flux  with  the  iono- 
spherically  induced  frequency  deviations. 


i\.ey  Words:  ionosphere,  solar  flare,  solar  X-ray  emission,  sudden  fre¬ 
quency  deviation  (SFD),  sudden  ionospheric  disturbance, 
solar  radio  emission 


A  COMPARISON  OF  SUDDEN  IONOSPHERIC  FREQUENCY  DEVIATIONS 
WITH  SOLAR  X-RAY  AND  CENTIMETER -WAVE  EMISSION  DURING  OCTOBER  1963 

*>y 

Donald  M.  Baker 

1.  Introduction 

This  report  presents  a  study  of  the  relationship  between  the  radio 
and  X-ray  emissions  of  solar  flare-  unu  the  ionospheric  disturbances 
called  sudden  frequency  deviations  (SFD).  Since  the  X-ray  bursts  which 
accompany  many  solar  flares  are  absorbed  in  the  earth's  atmosphere,  we  car 
use  ground-based  radio  techniques  to  detect  these  X-ray  bursts.  Hcwever, 
in  the  past  few  years  solar  X-ray  emission  has  been  monitored  directly  by 
several  artificial  satellites  (Bowen  et  al.,  Ih64;  Chubb,  Friedman,  and 
Kreplin,  1964;  Conner  et  al.,  1964;  Friedman,  1964;  Pounds,  1965 ) -  These  satel 
lite  observations  offer  us  an  opportunity  for  studying  the  relationship  between 
the  ionizing  radiation  and  the  resulting  disturbances  of  ionospheric  radio 
propagation,  such  as  sudden  phase  anomalies  (SPA),  short  wove  1  >deouts  (SWF), 
and  sudden  frequency  deviations  (SFD).  Such  studies  are  needed  if  ground- 
based  radio  techniques  are  to  be  used  as  tools  to  investigate  the  physics 
of  the  ionospheric  response  to  solar  ionizing  radiations  or  if  these 
techniques  are  to  be  used  to  study  the  nature  of  the  radiation  bursts 
themselves . 

There  is  evidence  that  solar  radio  bursts  at  centimeter  wavelengths 
and  X-ray  bursts  are  very  closely  related  (Kundu,  196h).  This  makes  a 
comparison  of  the  ionospheric  effects  of  solar  flares  with  both  the  X-ray 
and  associated  centimeter -wave  radio  emissions  desirable. 


In  this  report  we  will  restrict  ourselves  to  a  study  of  the  re. at  ion- 
ship  between  sudden  frequency  deviation  and  the  associated  X-ray  and  centi¬ 
meter  radio  bursts.  In  section  2  we  develop  a  simple  theory  useful  in  the 
interpretation  of  a  sudden  frequency  deviation  and  use  this  theory  to 
calculate  the  shapes  of  the  frequency  deviations  which  would  be  produced 
by  various  time  pr: files  for  the  ionizing  radiation  burst.  In  section  j 
we  compare  the  sudden  frequency  deviations,  centimeter-wave  radio  emis¬ 
sions,  and  the  X-ray  bursts  observed  by  the  Vela  satellites  for  the  six 
s.  Lur  flares  listed  in  table  1.  A  description  of  the  data  used  is  given 
be i ow . 

Table  1 

Optica  1  (:lor )  Observations  of  3,  lar  F.ares  Discussed  in  Se  -t :  n  * 

(E  -  before,  D  -  later  than) 

Time  (IT) 

Sate  beginning  Maximum  End  Lmportan 

Oct  .  lb.  2-uo  .  \*o  1  + 

1  <  .o'  ■  10  V  D 


1.2.  Centimeter  Radio  Data 


The  radio  data  consist  of  fixed  frequency  observations  at  centimeter 
wavelengths  made  at  Ottawa,  Ontario  (28 00  Me/s)  (Covington  and  Harvey, 

1958;  Harvey,  1964)  and  Toyokawa,  Japan  (200,  3750,  and  9400  Mc/s) 

(Tanaka  and  Kakinuma,  1958).  The  events  used  in  this  report  have  been 
replotted  from  copies  or  tracings  of  the  original  records. 

1.3.  Ionospheric  Disturbance  Data 

The  ionospheric  disturbance  data  consist  of  variations  in  the  received 
frequency  of  high- frequent /  ion  spheri  all’  propagated  radio  signals.  These 
data  are  obtained  by  the  Doppler  technique  described  by  Watts  and  Davies 
(i960)  which  detects  rapid  changes  in  the  electron  content  of  the  iono¬ 
sphere  up  to  the  level  of  reflection. 

The  frequency  of  an  ionospherically  propagated  HF  radio  signal 
usually  shows  small  variations  (of  the  order  of  a  few  tenths  of  a  cycle 
per  second)  about  the  transmitted  frequency.  However,  during  some  natural 
phenomena,  such  as  solar  fiares  and  geomagnetic  sudden  commencements,  the 
received  frequency  often  shows  large  (up  to  tens  of  cycles  per  second) 
and  quite  distinctive  variations  from  the  transmitted  frequency.  These 
flare-related  frequency  deviations  have  been  given  the  name  sudden  fre¬ 
quency  deviations  (SFD)  (Chan  and  Villard,  1963 ) *  btuaies  by  Donnelly 
(1966);  Agy,  3aker,  and  Jones  (1965);  Davies,  Watts,  and  Zacharisen  (1962); 
and  Kanel 1  kcs,  Chan,  and  Villard  (1962)  indicate  that  they  are  caused  by 
an  increase  of  ionization  in  the  E  and/or  F  regions  of  the  ionosphere. 
Details  of  the  technique  and  the  interpretation  of  the  records  can  be 
found  in  Dav'es  and  Baker  (1966)  and  Donnelly  (i960)  and  the  references 
contained  therein. 
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2.  Theoretical  Discussion 
2.1.  Basic  Equations  and  Model  Used 
If  a  radio  wave  of  carrier  frequency  f  propagates  through  a  changing 
ionosphere,  it  will  suffer  a  change  in  frequency,  or  Doppler  shift,  Af 
given  by 


Af  =  - 


f 

c 


dP 

dt 


where  c  is  the  speed  of  light  in  vacuum.  The  phase  path  of  propagation, 
P,  is  given  by 

P  =  /  Hds  ,  (2) 

path 


where,  neglecting  the  effects  of  the  geomagnetic  field  and  collisions, 

A  =  V  1  '  p  (3) 

is  the  refractive  index.  Here  k  =  8  x  10"6  (?)•  cm3  ,  the  electron  den¬ 
sity  N  is  in  cm'3,  and  the  wave  frequency  f  is  in  Mc/s. 

Using  (2)  and  (3),  Agy,  Baker,  and  Jones  (i960)  have  shown  that  (l) 
can  be  written 


Af  * 


k 

2fc 


ds  . 


(k) 


If  the  changes  causing  Af  are  confined  to  a  non-deviating  region 
where  p  «»  1  (i.e.,  a  region  below  the  height  of  reflection) 
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reflection 


grouna 


kn  dNr 
fc  at 

where 


(5) 


dNr 

dt 


reflection 


/ 


dN 

dt 


ds 


ground 


(6) 


is  the  rate  of  change  of  the  total  electron  content  and  n  is  the  number  of 
ionospheric  reflections.  Therefore,  when  ionospheric  changes  are  confined 
to  a  region  below  the  reflection  height  of  the  radio  wave  being  monitored, 
the  Doppler  shift  observed  will  be  directly  proportional  to  the  total  time 
rate  of  change  of  electron  content  along  the  propagation  path.  Ii  the 
ionospheric  changes  extend  all  the  way  to  the  height  of  reflection  and  if 
there  is  no  major  change  in  the  propagation  path  (such  as  a  change  from 
F-layer  to  E-layer  reflection),  we  would  expect  from  (4)  that  the  shape  of 
Af(t)  would  still  reflect  the  major  time  variations  of  the  total  electron 
content  along  the  path;  however,  the  magnitude  of  Af(t)  will  now  be  influ¬ 
enced  by  the  variation  of  p  along  the  path.  Henceforth,  we  will  consider 
only  the  non-deviative  model. 

The  time  rate  of  change  of  the  electron  density  in  a  region  in  which 
the  electron  loss  processes  cud  be  described  by  an  effective  recombination 
coefficient  b  is  given  by 


dN 


dt  ^ 


air  ,  (7) 

where  q  is  the  electron  production  rate.  Assume  that  the  electron  pro- 
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duction  and  loss  rates  are  in  equilibrium  before  a  solar  flare.  Then,  if 
the  production  rate  and  electron  density  are  increased  by  amounts  Aq(t) 
and  AN(t)  during  a  flare  (i.e.,  q(t)  =  q0  +  Aq(t),  N(t)  =  N0  +  AN(t)  where 
the  zero  subscripts  denote  the  equilibrium  values),  (7)  can  be  written  as 


or 


where  t  is  defined  as 


1 

2oNu  * 


(8) 


(9) 


(10) 


Integration  of  (9)  over  the  same  integration  path  used  to  obtain  (5), 
assuming  that  the  region  affected  by  the  flare  can  be  characterized  by 
constant  a  and  NQ,  yields 

mill  .  lq,(t)  -  "!lM  (1  ,  Ml gl\  ,  (u) 

where  the  time  dependence  is  now  indicated  explicitly.  This  equation, 
when  combined  with  ([>),  allows  us  to  synthesize  the  Af  to  be  expected 
from  a  known  or  assumed  Aqr  (t),  or  to  deduce  the  Aqj  (t)  responsible  for  an 
observed  Af(t).  Moreover,  since  the  frequency  deviations  caused  by  flares 
rarely  Inst  more  than  a  few  minutes,  the  solar  zenith  angle  is  essentially 
constant  during  such  an  event.,  and  the  electron  production  rate  in  a  given 
region  should  closely  follow  the  flux  of  the  ionizing  radiation  enhancement. 
Hence,  if  the  major  flare  effects  detected  by  the  Doppler  technique  do 
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indeed  take  place  below  the  height  of  reflection,  the  simple  model  devel¬ 


oped  here  should  permit  us  to  determine  the  time  profile  of  the  burst  of 
ionizing  radiation  causing  an  ionospheric  disturbance. 

When  AN  «  2NQ,  (ll)  can  be  solved  analytically  for  dANr/dt.  'This  is 
most  easily  done  by  differentiating  (ll)  with  respect  to  time  and  solving 
the  resulting  second  order  equation  for  A N(t),  where  we  now  let  the  dot 
stand  for  the  time  derivative.  The  result  gives: 


t 

ANr  (t)  =  e  T 


Aqt  (t  )dt 


(12) 


where  C  is  a  constant  of  integration.  This  last  expression  can  be  stated 
in  terms  of  Aqt  (t)  instead  of  A4t  (t)  by  integrating  by  parts: 


ANt  (t)  =  AqT  (t)  - 


j AqT  (t ) 


dt 


(13) 


If  Aqr  (t)  and  the  appropriate  value  of  t  are  known.  (13)  enables  ANt  (t), 
and  hence  Af(t),  to  be  calculated.  Under  these  conditions  (i.e.  AN  «  2Nq), 
t  can  be  considered  as  an  effective  relaxation  time  of  the  ionosphere. 

Alternatively,  the  AqT  (t)  responsible  for  an  observed  Af(t)  can  be 
obtained  by  substituting  (‘ )  into  (ll)  to  obtain  the  following  expression 
for  AqT  (t)  in  terms  of  Af(t): 

,(U0 

where  x  i  s  a  vur  i  ab  1  c*  of  i ntograt  i on . 
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AqT  (t  ) 


fe 

k 


Af(t) 


i.  _  t 

7 /  M'(x)ox  i‘|i/‘r(-l'K 

o  L  "  ° 


$ 


Y 


Examples  of  both  of  these  approaches  will  be  given  in  section  3,  where 
the  frequency  disturbances  are  compared  with  the  X-ray  bursts  in  the  0.5- 

o 

to  10-A  ranges  observed  by  the  Vela  satellites. 

2.2  Model  Calculations 

Five  simple  models  of  the  enhancement  of  the  electron  production  rate, 
Aqr  (t),  are  shown  in  figures  l(a),  l(c),  and  1(e)  and  in  figures  2(a)  and 
2(c).  The  frequency  deviations  caused  by  these  production  rate  enhance¬ 
ments  were  calculated  from  (13)  using  assumed  values  of  the  effective 
relaxation  time,  t,  of  0.5,  1,  5,  and  10  min.  Tnree  simple  models  of 
Aqr  (t)  are  shown  in  figure  1.  An  impulsive  rise  to  a  peak  followed  by  a 
slower  decay  to  the  undistrubed  level  (fig.  la)  causes  a  positi/e  fre¬ 
quency  deviation  followed  by  a  negative  shift  and  a  gradual  recovery  to 
the  undisturbed  frequency  (fig.  lb).  A  rapid  rise  of  electron  production 
rate  to  a  new  constant  level  (fig.  lc)  causes  a  rapid  frequency  increase 
and  a  gradual  recovery  to  the  undisturbed  frequency  with  no  negative 
phase  (fig.  Id).  A  gradual  rise  and  fall  of  the  electron  production  rate 
(fig.  le)  results  in  a  small  and  gradual  increase  in  frequency  followed 
by  a  small  negative  deviation  and  a  gradual  recovery  (fig.  If).  Fre¬ 
quency  deviations  similar  to  those  in  figures  i(b)  and  l(d)  have  been 
observed  frequently  during  solar  flares.  Effects  similar  to  those  of 
figure  l(f)  have  also  been  observed;  hcvever,  such  small,  gradual  fre¬ 
quency  variations  are  hard  to  distinguish  from  the  normal  background  vari- 
3ti  ons  which  are  usually  present,  and  up  to  now  they  have  not  been  con¬ 
sidered  as  flare-related  events.  In  the  past  the  rapidity  of  a  frequency 
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deviation  has  been  one  of  the  criteria  used  to  recognize  a  sudden  frequency 


deviation;  the  present  study  suggests  that  such  a  criterion  may  result  in 
many  flare-related  events  being  missed.  However,  the  Doppler  technique 
is  not  suitable  for  quantitative  study  of  these  gradual  events. 

Slightly  more  complex  electron  production  rate  enhancements  are 
shown  in  figures  2(a)  and  2(c).  Figure  2(a)  shows  a  Aq  with  a  double  rise, 
the  two  increases  being  separated  by  a  short  period  during  which  the  pro¬ 
duction  rate  remains  at  a  constant  but  enhanced  level.  For  the  shorter 
relaxation  times  such  a  model  leads  to  a  frequency  deviation  having  two 
distinct  peaks  (fig.  2b),  which  become  less  pronounced  as  the  relaxation 
time  increases.  An  impulsive  burst  superimposed  upon  the  initial  part  of 
a  gradual  rise  and  fall  is  shown  in  figure  2(c).  This  production  rate 
model  causes  an  initial  frequency  deviation  similar  to  that  caused  by  a 
simple  impulsive  burst  (fig.  la  and  lb);  however,  for  the  shorter  relaxa¬ 
tion  times  the  graduaL  increase  in  the  production  rate  following  the  im¬ 
pulsive  burst  shortens  the  length  of  the  negative  deviation  and  causes  a 
sec.  rid  p  s Stive  frequency  shift  similar  to  that  caused  by  a  simple  gradual 
v ; s-  and  fa. 1  of  the  produce i  n  rate  (fig.  ie  und  If).  This  "overshoot" 
is  full  wet  by  a  small  negative  deviation  and  a  gradual  recovery.  For  the 
longer  ro'.axati  n  times  the  grades .  rise  of  the  production  rate  prevents 
he  frep.eney  devia*.  i;  n  fr  m.  reaching  the  initial  negative  phase. 

Many  .  f  *h"  f la-e-re  fated  frequency  disturbances  observed  have  more 
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smaller  than  .  r  mparable  t_  the  time  between  rapia  flux  changes  (about 


1  min  in  this  example)  the  burst  of  ionizing  radiation  need  not  have  dis¬ 
tinct  peaks  but  merely  needs  to  be  made  up  of  rapid  flux  increases  sepa¬ 
rated  by  plateaus  in  order  to  cause  distinct  peaks  in  the  frequency  deyia- 
tion.  Frequency  disturbances  made  up  of  a  main  deviation  followed  by  an 
overshoot  such  as  those  of  figure  2(d)  have  also  been  observed.  Clearly, 
such  events  could  be  caused  by  an  impulsive  burst  of  ionizing  radiation 
superimposed  upon  a  gradual  rise  and  fall  similar  to  the  production  rate 
model  of  figure  2(c). 

Comparison  of  the  frequency  deviations  for  the  different  relaxation 
times  and  the  appropriate  enhancements  of  the  electron  production  rate  in 
figures  1  and  2  shows  (see  (9))  that  the  frequency  variation  follows  the 
enhancement  of  electron  production  rate  more  closely  as  the  effective 
relaxation  time  becomes  larger.  The  frequency  deviation  can  depart  sig¬ 
nificantly  from  the  shape  of  the  enhancement  of  the  production  rate  for 
small  relaxation  times.  Also,  a  negative  frequency  deviation  is  to  be 
expected  only  when  the  production  rate  (i.e.,  the  ionizing  flux)  decreases 
with  time. 


i.  Specific  Events 

In  this  section  we  compare  the  X-ray  bursts,  radio  bursts,  and  fre¬ 
quency  deviations  observed  during  the  flares  listed  in  table  1.  The  times 
of  onset,  durations,  and  general  behavior  of  the  three  phenomena  are  inves¬ 
tigated  first.  The  shape  of  the  frequency  deviation  to  be  expected  from 
•he  0.5-  to  10-X  X-ray  burst  of  October  22  is  then  calculated.  Finally, 
the  time  profiles  of  the  enhancements  of  electron  production  rate  are 
deduced  and  compared  with  thr  X-ray  observations. 


3.1.  Comparison  of  Observed  X-ray  Bursts,  Gclar  Radio  Bursts, 
and  Frequency  Deviations 

The  time  variations  of  the  X-ray  flux  in  the  0.5-  to  4-1  and  0.5-  to 
10-1  ranges,  the  solar  radio  flux  at  fixed  frequencies  in  the  centimeter 
wavelength  region,  and  the  frequency  deviation  of  WWV  (or  WWVH)  as  received 
at  Boulder  are  shown  in  figures  3  through  8.  The  radio  flux  is  given  in 
flux  units  (l  flux  unit  =  10  22  ,  W  m  2  cps'1 ),  the  frequency  deviation 
in  cycles  per  second,  and  the  X-ray  flux  deposited  in  the  detector  in 
ergs  I’m"9  sec"'  .  Note  that  the  X-ray  flux  scale  is  linear.  The  0.5-  to 
10-Jl  detectors  became  saturated  during  all  but  the  2240  event  uf  October  22 
the  0.5-  to  4-4  detectors  were  saturated  only  during  the  October  18  event. 
The  4- min  gaps  in  the  Doppler  data  which  occur  at  4*;  min  past  the  hour 
are  due  to  the  hcu.  ly  WWV  transmission  break.  The  lowest  level  shown  for 
the  X-ray  flux,  about  1  x  10-3  ergs  cm-3  sec-1 ,  represents  the  detector 
threshold  and  is  not  necessarily  the  background  level.  The  radio  and 
D.  ppier  data  have  been  scaled  and  replotted  from  copies  of  the  original 
re  u  r  is;  the  t  iming  is  judged  to  be  accurate  to  within  a  minute  or  better. 
The  ’  iming  .  f  the  plots  .  f  the  X-ray  data  may  be  off  as  much  as  5  mi*,. 

Cum; aris-'n  f  the  observed  0.5-  to  I0-A  X-ray,  centimeter  radio,  and 
fr-*q  variati.n  data  ?  r  the  six  events  leads  to  the  following  general 

a*  ns: 


T.  wi'oin  *  h«*  timing  accuracy  available,  the  unset  times 
:  the  Imjulsive  radi  bursts  and  the  frequency  deviations 
are  in  agr-em.en* . 

The  X-ray  flux  usually  increases  rapid,  y  iu*  ing  the  impul- 
ve  rauij  bursts  and  the  frequency  deviations.  The 


deviations. 


Doppler  technique  tends  to  "see"  the  initial  fast  rise  of 
the  X-ray  flux. 


3.  In  the  case  of  four  of  the  five  events  for  which  radio 
data  are  available  the  duration  of  the  frequency  devia¬ 
tion  is  comparable  to  that  of  the  impulsive  radio  burst; 
for  the  fifth  event  (fig.  5),  the  impulsive  radio  burst 
lasts  longer  than  the  frequency  deviation. 

4.  For  some  events  (figs.  3  and  6)  the  duration  of  the 
X-ray  burst  is  of  the  same  order  as  that  of  the  fre¬ 
quency  deviation;  during  other  events  (figs.  4,5,7, 
and  8)  the  X-ray  flux  remains  at  an  enhanced  level 

long  after-  the  frequency  disturbance  has  ended. 

5.  No  significant  frequency  deviations  are  observed  during 

the  (apparently)  rapid  decay  of  a  long  enduring  X-ray 
burst.  However,  the  exact  beginning  time  and  rate  of 
decay  of  the  recovery  phase  of  such  3  burst  c..nnct  be 
determined  due  to  the  sc*urat l on  ef  the  detectors. 

6.  The  radio  bursts  and  frequency  deviations  often  show 
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j.2.  Calculation  of  the  Frequency  Deviation  Expected  from  the 
X-ray  Burst  Observed  at  2240  U.T.  on  22  October  19^3 


Equation  (13)  of  section  2  permits  us  to  calculate  the  frequency  de¬ 
viation  which  will  be  caused  by  a  given  enhancement  of  the  electron  pro¬ 
duction  rate.  In  deriving  (13)  we  assume  that  the  changes  responsible  for 
the  frequency  deviation  are  confined  to  a  limited  portion  of  the  iono¬ 
sphere  which  is  nrn-devia uive  for  the  radio  frequen  ces  used.  This  effec¬ 
tively  limits  the  changes  to  the  E  region  and  lower  p  rt  of  the  F  region 
for  observations  made  at  10  and  15  M'  s.  We  also  assune  that  the  electron 
loss  processes  in  *his  regi.n  can  be  repr-.sented  by  an  effective  relaxa¬ 
tion  time,  t ,  which  is  essentially  constant  with  height  and  for  times  of 
the  order  of  the  duration  f  the  frequency  deviation  of  interest.  For  a 
region  where  re . 'cabinet  1. n  Css  pr  cesses  predominate,  this  assumption 
requires  that  the  increase  in  e.ectrcn  density  caused  by  the  flare  be 

sme  1 1  c  mpareu  with  the  ambient  election  density.  Donnelly  (:96b)  has 

* 

si.  vn  that  X-radiati  n  a*  wavelengths  ar  uni  if  A  : reduces  ionization 
m  s*  iy  in  *h~  S'-  *  •  if- Km  regi.n.  This  leve.  is  we  1 1  be.cv  the  height 

f  ref.ee*  i.  n  (i.e.,  is  n  n-devia*  ive )  f.  r  F-.sy-r  ;r  oogati  n  of  10  to 
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quency  deviation  tc  be  expected  from  a  given  X-ray  burst,  we  will  assume 
that  the  enhancement,  of  the  electron  production  rate  has  the  same  time 
profile  as  the  X-ray  burst. 

The  results  of  applying  this  technique  to  the  0.5-  to  10-A  X-ray 
event  'braved  by  the  Vela  satellites  at  22h0  U.T.  on  22  October  19&3  are 
shown  in  figure  9*  The  "stair-step"  profile  of  th-  burst  shown  in  figure 
6(c)  has  been  smoothed  as  shown  in  figure  9(8)*  This  profile  was  then 
used  for  Aqr(t)  in  (13),  and  the  frequency  deviation  was  calculated  for 
effective  relaxation  times  of  0.5,  and  5  min.  The  results  are  shown 

by  the  dotted  curves  of  figures  '  9(c),  and  9(d);  the  solid  curve  shows 

the  observed  frequency  shift  for  comparison.  The  times  of  the  calculated 
and  observed  peaks  have  been  aligned  and  the  calculated  peaks  have  been 
normalized  to  agree  with  the  maximum  observed  deviation.  There  is  some 
evidence  that  the  X-ray  peak  occurred  about  1  min  after  the  peak  fre¬ 
quency  deviation  (see  fig.  12);  however,  this  difference  is  well  within 
the  accuracy  with  whi  h  the  rec.rds  are  timed  and  for  the  purposes  of  this 
calculation  the  teaks  were  assumed  to  occur  at  the  same  time. 

From  figures  ‘■'(c)  and  9(d)  we  see  that  the  gross  features  of  the  ob¬ 

served  and  calculated  frequency  devir*  i ons  are  in  good  agreement  for  re- 
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satellites  cow  1.1  have  produced  the  .  fc.««rved  frequency  deviation.  We  have 
said  nothing  ab  u*.  the  magnitudes.  In  a  more  detailed  study  f  this  event. 
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Donnelly  (19 66)  has  found  that  the  magnitude  of  the  flux  enhancement 
in  the  0*5-  to  10-A  range  is  not  sufficient  to  produce  the  observed  fre¬ 
quency  deviation.  He  concludes  that  the  frequency  deviation  must  be 
at  least  partly  a  result  of  enhancements  at  longer  wavelengths. 

Since  the  flux  exceeded  the  saturation  level  during  a  major  portion  of 
the  other  events  observed  at  Boulder,  it  is  not  worthwhile  to  reconstruct 
the  flux  profile  and  calculate  the  frequency  deviations  which  would  be 
produced. 


3.3.  Synthesis  of  Profiles  of  the  Enhancement 
of  Electron  Production  Rate 

In  this  section  we  take  the  observed  frequency  deviations  as  the  input 
data  and  use  (l4)  to  calculate  th-  enhancement  of  the  electron  production 
rate.  The  resulting  profiles  for  effective  relaxation  times  of  0.5,  1,  and 
10  min,  shown  in  part  (b)  of  figures  10  through  14,  are  compared  with  the 
observed  0.5-  to  10-A  X-ray  flux  shown  in  part  (c)  of  these  figures.*  The 
time  scale  of  these  figures  is  much  expanded  over  that  used  in  figures  3 
through  8,  and  only  the  initial  portion  of  each  X-ray  event  is  shown.  After 
the  initial  peak  frequency  deviation  and  the  initial  part  of  the  negative 
shift,  the  flare-related  frequency  deviation  is  usually  lost  in  the  back¬ 
ground  variation  of  the  frequency;  in  making  the  calculations  for  figures 
10  through  14,  only  that  portion  of  a  frequency  deviation  which  could  be 
reliably  attributed  to  the  flare  was  used. 

*In  calculating  the  enhancement  of  the  electron  production  rate  we  have 
used  a  value  of  1011  m“3  for  N  , 
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The  curves  of  figures  10  through  14  show  the  shape  of  the  electron 
production  rate  enhancement  which  would  be  needed  to  produce  the  observed 
frequency  deviation  for  the  simple  model  we  have  adopted.  Comparison  of 
the  production  rate  enhancements  with  the  X-ray  fluxes  must  of  necessity 
be  qualitative  due  to  the  nature  of  the  data  available;  our  knowledge  is 
limited  by  the  saturation  of  the  X-ray  detectors,  the  lack  of  fine  struc¬ 
ture  in  the  X-ray  observations  and  the  limited  portion  of  the  production 
rate  profile  which  we  can  reconstruct  from  the  frequency  variation  data. 
The  most  we  can  hope  to  conclude  is  whether  or  not  the  calculated  produc¬ 
tion  rate  enhancement  is  compatible  with  the  X-ray  observations. 

Smith,  Ac car do,  Weeks,  and  McKinnon  (1965)  have  concluded  from  recent 
eclipse  observations  that  the  effective  daytime  relaxation  time  in  the  E 
region  is  about  1  min.  In  addition,  baker  and  Davies  (1966)  have 
presented  evidence  that  suggests  that  the  effective  daytime  relaxatxon 
time  in  the  region  where  the  bulk  of  the  ionization  responsible  for  flare- 
related  frequency  deviations  is  released  may  be  less  than  1  min. 

Therefore,  the  most  appropriate  curves  for  the  enhancement  of  electron 
production  rate  for  the  model  we  have  adopted  are  probably  those  for  an 
effective  relaxation  time  of  about  1  min. 

Examination  of  figures  10  and  11  reveals  that  the  synthesized  pro¬ 
duction  rate  enhancement  curves  and  the  initial  part  of  the  X-ray  burst 
are  compatible  for  these  two  events.  From  a  cursory  inspection  it  may 
appear  that,  the  production  rate  curve  is  decaying  faster  than  the  X-ray 

flux,  but  closer  inspection  reveals  that  at  the  end  of  the  records 

shown  in  the  figures,  the  production  rates  are  still  above  the  levels  they 


had  attained  at  the  time  the  X-ray  detectors  became  saturated.  The  pro¬ 
duction  rate  curves  for  these  two  events  also  show  that  the  fine  structure 
in  the  ionizing  radiation  needed  tc  produce  the  fine  structure  of  the  fre¬ 
quency  deviations  is  indeed  small,  especially  if  the  effective  relaxation 
time  is  short,  and  such  small  relative  flux  variations  may  prove  to  be 
hard  to  detect  with  satellite-borne  instrumentation. 

Figure  12  shows  the  corves  for  the  2240  October  22  event  which  was 
discussed  in  section  3-2.  The  time  disagreement  between  the  X-ray  burst 
and  the  frequency  deviation  which  was  mentioned  previously  is  shown  in 
this  figure.  Here  we  will  adopt  the  viewpoint  that  this  discrepancy  is 
within  the  accuracy  to  which  the  data  can  be  timed  and  will  assume  that 
the  two  events  are  time  coincident.  Doing  so,  we  see  that  the  production 
rate  curves  for  the  shorter  relaxation  times  are  not  compatible  with  the 
X-ray  observations;  indeed,  the  production  curves  decay  much  more  slowly 
than  the  X-ray  flux.  The  production  rate  curve  for  a  relaxation  time  of 
10  min  is  in  fairly  good  qualitative  agreement  with  the  X-ray  flux 
curve;  this  we  would  expect  from  the  results  of  section  3-2  where  we 
found  that  a  relaxation  time  greater  than  1  min  was  needed  to  make 
the  frequency  deviation  calculated  from  the  X-ray  flux  agree  with  the 
observed  event. 

The  initial  part  of  the  event  of  October  26  is  shown  in  figure  13. 

Only  the  two  prominent  frequency  deviations  at  1840  and  1902  were  consid¬ 
ered  in  calculating  the  enhancement  of  the  electron  production  rate;  the 
data  between  these  two  deviations  were  smoothed  as  indicated  by  the  f 

dashed  lines.  From  the  portion  of  the  event  shown,  the  production  rate  | 
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enhancement  and  the  X-ray  burst  appear  to  be  compatible.  However,  the 
X-ray  detectors  had  not  yet  reached  saturation  at  1908,  and  an  abrupt 
increase  to  the  saturation  level  shortly  after  1908  caused  no  measurable 
frequency  deviation.  This  fact  forces  us  to  question  whether  the  0.5"  to 

o 

10-A  X-ruy  flux  was  responsible  for  the  frequency  deviations  before 
1908. 

The  electron  production  rate  enhancement  calculated  for  October  28 
(fig.  14)  does  not  agree  in  time  with  the  X-ray  flux  enhancement.  The 
production  rate  has  peaked  and  is  beginning  to  deca^  before  the  X-ray 
detectors  become  saturated,  A  shift  of  about  4  min  would  be  necessary 
to  make  the  production  rate  curve  compatible  with  the  X-ray  burst; 
again  such  a  shift  is  within  the  timing  accuracy  available. 

Prom  figures  10  through  l4  we  see  that  in  some  cases  the  X-ray  flux 
in  the  0.5-  to  10-A  range  could  have  caused  an  enhancement  of  the  electron 
production  rate  of  the  shape  needed  to  give  the  observed  frequency  varia¬ 
tion,  and  that  in  other  cases  it  could  not.  It  may  be  significant  that 
for  the  one  event  during  which  the  X-ray  detectors  were  not  saturated 
(fig.  12),  the  production  rate  enhancement  calculated  for  the  shorter 
relaxation  times,  which  we  would  expect  to  be  operative  in  the  E  region 
where  our  model  is  applicable,  differs  appreciably  from  the  X-ray  burst. 
As  pointed  out  by  Donnelly  (1966) ,  th«  o  ^  10-1  flux  may  have  contri¬ 

buted  to  the  post -peak  hump  in  the  frequency  variation  while  the  main 
deviation  was  caused  by  enhancements  at  wavelengths  greater  than  10  A. 


4.  Discussion 


We  cannot  establish  by  a  comparison  of  the  time  profiles  of  the  data 
available  whether  or  not  X-ray  bursts  in  the  0.5-  to  10-1  range  were  suf¬ 
ficient  to  cause  the  frequency  deviations  observed  during  the  solar 
flares  of  October  1963.  To  compare  the  frequency  deviations  and  X-ray 
data,  we  have  assumed  that  the  spectral  composition  of  the  X-ray  burst 
did  not  change  during  the  burst.  Some  of  the  disagreement  we  have,  found 
may  arise  because  the  spectral  characteristics  of  an  X-ray  burst  do  change 
rapidly  with  time.  In  any  case,  we  could  not  definitely  conclude  that  the 
frequency  variations  were  or  were  not  caused  by  the  0.5-  to  10-1  X-rays 
from  a  study  of  time  profiles  alone  without  knowledge  of  the  flux  varia¬ 
tions  at  longer  wavelengths.  Donnelly  (19 66)  has  concluded  that  the  fre¬ 
quency  deviations  must  be  at  least  partly  a  result  of  enhancements  at 
wavelengths  greater  than  10  A. 

Such  a  simple  treatment  of  the  data  does  reveal  some  worthwhile  ob¬ 
servations,  however.  We  see  that  better  timing  accuracy  is  required  for 
detailed  comparison  of  the  various  data;  such  timing  accuracy  is  now 
available  on  the  Doppler  data.  Attempts  should  be  made  to  obtain  more 

detailed  time  profiles  and  better  spectral  resolution  of  solar  X-ray 
bursts.  Although  the  calculated  enhancements  of  the  electron  production 
rate  indicate  that  the  fine  structure  necessary  to  produce  the  detai 1 
observed  in  the  frequency  deviations  may  be  -nly  a  very  small  modulation 
of  the  main  enhancement,  the  elose  correspondence  which  often  exists 
between  the  frequency  variations  and  the  impulsive  centimeter  radio  bursts 
suggests  that  similar  impu’sivo  bursts  may  "xist  in  the  flux  of  ionizing 
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radiation.  Better  spectral  resolution  and  broader  spectral  coverage  are 
needed  to  help  us  determine  what  wavelength  regions  produce  the  ionization 


enhancements  responsible  for  the  frequency  variations  and  the  height  of 
these  enhancements. 
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The  solar  radio  burst  (a),  frequency  variation  (b),  and  solar 
burst  (c)  for  the  event  of  l8  October  1963  (2041?  U.T.).  Ail  , 
are  linear. 


<Q>  SOLAR  RADIO  BURST,  2600  MC/S ,  OTTAWA 


1630  UT 


!  solar  radio  burst  (a),  frequency  variation  (b),  and  solar  X-ray  burst  (c)  for  the  event  of 
October  1963  (1330  U.T.).  All  scales  are  linear. 
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(a)  SOLAR  RADiO  BURSTS,  TOYOKAWA,  JAPAN 


2000  MC/S 


Fig.  6  The  solar  radio  burst  (a),  frequency  variation  (b),  and  solar 
X-ray  burst  (c)  for  the  event  of  22  October  1963  (2240  U.T.). 
All  scales  are  linear. 
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(a!  SOLAR  RADIO  BURST ,  2800MC/S,  OTTAWA 
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The  frequency  deviation  and  solar  X-ray  burst  for  the  event  of  28  October  1963  (0140 
scales  are  linear. 


Fig.  9  The  0.5-  tc  10-A  X-ray  flux  (a)  and  the  synthesized  (dotted)  and 
observed  (solid  line)  frequency  variations  (b),  (c),  and  (d)  for 
the  event  of  22  October  1963  (2240  U.T.).  The  peaks  of  the  syn¬ 
thesized  and  observed  frequency  variations  have  been  normalized 
and  aligned. 
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FREQUENCY  VARIATION,  1 5MC/S,  W*V  TO  BOULDER 


Fi £ .  12  The  observed  frequency  variation  (a)  and  X-ray  burst  (c)  and  the 
•aiculatej  enhancement  of  the  electron  production  rate  for  the 
event  of  ^  October  196  j  (22^0  U.T.).  All  scales  are  linear. 

3b 


(o)  FREQUENCY  VARIATION,  i5Mc/S,  WWVH  TO  BOULDER 
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